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Abstract 

We have proposed a model of roton cluster excitations in liquid helium II based on a Schrodinger- 
type equation with a self-consistent confining potential. We have derived an equation for the 
number of atoms in roton excitations, which can be treated as quantum 3D solitons, depending 
on vibrational quantum numbers. It is shown that the smallest roton cluster is in the symmetric 
vibrational quantum state and consists of 13 helium atoms. We have also used a modified Born 
approximation to calculate the s-scattering length for helium atoms. This allows us to calculate all 
parameters of Landau's roton excitation spectrum, in agreement to high accuracy with experimental 
measurements from neutron scattering. 

PACS numbers: 67.40.Db, 03.75.Fi, 05.30.Jp 
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1 INTRODUCTION 



Liquid helium is a unique system of strongly interacting atoms which becomes a quantum 
liquid (helium II) when the temperature is lower then some critical temperature T c . This 
critical temperature T c can be well estimated by two parameters: the de Broglie thermal 
wavelength A^ = 2irh(3'mkBT)~ 1 / 2 and the s-scattering length ao = 2.2 A of the helium 
atoms. Thus we may suppose that in liquid helium, a phase transition occurs to a quantum 
liquid when the de Broglie thermal wavelength overlaps two nearest neighbors of helium 
atoms: \d — 2c?, where d = 2ao is the effective diameter of helium atoms. This equation 
yields the temperature T c = 2.05 K, which is very close to the experimental value T c = 2.18 K 
of the phase A-transition in liquid helium. 

The study of liquid helium II, and quantum liquids more generally, has been an active 
area of experimental and theoretical research ever since Kapitza's [l| discovery of the su- 
perfluidity in helium II. Based on some experimental facts, Landau Q] assumed that in 
liquid helium II there exist two types of elementary excitations: phonons and rotons. He 
also phenomenologically introduced the roton energy spectrum (see Eq.(£Q)) to explain the 
thermodynamic behavior of liquid helium II. Bogoliubov j^] first found theoretically the 
spectrum of quasi-particles in a dilute Bose gas. This theory can be used for a Bose gas 



when the temperature is close to zero, the particle density n is low (\J a^n <C 1), and the 
wavenumbers of the particles are sufficiently small {ka^ <C 1). The relation between the 
energy spectrum of the elementary excitations of liquid 4 He and the structure factor of the 
liquid was found by Feynman \\ . Feynman also proposed a phenomenological model of the 
roton as a localized vortex ring with characteristic size of the order of the mean atomic 
distance in liquid helium II. The vortex rings were observed experimentally by Rayfield and 
Reif [sj]; however, there has been no experimental confirmation that the vortex ring and 
roton are the same excitations in liquid helium. 

We also note that a variety of many-body theories have been applied to the problem since 
1990 and a quantitative description of a roton has been obtained from variational integral 
equations, from a variational Monte Carlo method and by an imaginary time exact path 
integral. Thus, a correct model of roton quasi-particles should be based on "first principles" 
and lead to quantitative agreement with experimental data. Recently such a model was 
proposed in the paper p, where it was found that the roton is an excitation of a cluster 
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of 13 atoms, presumably having a central atom surrounded by a shell of 12 atoms situated 
at the vertices of a regular icosahedron. In that paper we considered the number of atoms 
iV c of the roton cluster as a fitting parameter that was found from the experimental data of 
Yarnell et al |7[. 

In fact, the theory developed in |(| does not need fitting to experimental data. To 
formulate this theory completely self-consistently we should find an additional equation for 
the chemical potential of roton excitations. In this paper, we calculate the vibrational energy 
spectrum of the roton cluster excitations with an arbitrary number of atoms in the cluster 



by the Schrodinger-type equation with a self-consistent confining potential derived in 
Using scattering theory we have formulated the modified Born approximation (MBA) and 
have found the s-scattering length a for helium atoms from the Lennard- Jones interatomic 
potential. We also derive in this paper the equation for the chemical potential \i of a roton 
cluster, which allows us to show that the smallest stable roton clusters in liquid helium II 
indeed contain 13 atoms. Moreover, we have found the numbers of atoms iV c for the roton 
clusters which are in a vibrational state \n\nin^ with arbitrary quantum numbers ni, ri2, n^. 
This self-consistent approach allows us to accomplish all calculations for the parameters 
of the roton cluter excitations in helium II and leads to a quantitative agreement with 
experimental data. In particular, the roton clusters with N c = 13 are in the symmetric 
vibrational state |001) and the energy spectrum of roton cluster excitations can be derived 
in the form [f| : 

£(p) = a+ (ezj*)!, (1) 

which coincides with Landau's empirical formula. 

2 SCHRODINGER EQUATION FOR ROTON EXCITATIONS 

The nonlinear Schrodinger equation (NLSE) for roton cluster excitations can be de- 
rived by a Hartree-Fock variation procedure for the trial wave function \I/_r(xi, ...,x n ;£) = 
C , 0ii(xi, t)... , 0a(x n , t) where n = N c . In this procedure we represent the interaction poten- 
tial by two parts: Ur(x.) and C/a(x), where £/#(x) is the repulsive part and Ua(x) is the 
attractive part. The repulsive part of this potential can be approximated by a hard sphere 
potential, 

£Mx-x') = ^5(x-x'), (2) 
m 
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and the attractive part (the confining many-body effective potential) is 

U A (x) = (m/2)j2"l4- (3) 

k=l 

The parameters uj\ are found by the stationarity conditions (dV/dak)a k =a k = (see below). 
This potential takes into account all forces between 4 He atoms of the cluster and those 
of the surroundings, and also long-range many-body attractive forces among atoms of the 
cluster This Hartree-Fock variation procedure directly leads to a NLSE : 

dip R ( h 2 2 m 3 2 2 Ana h 2 (N c - 1) 2 \ 
%h ^T = {-2^ V + 2 g^* + ^N c ^ ) (4) 

Here the normalisation condition for the wave function iPr(x., t) has the form f v |?/>#| 2 <i 3 x = 
N c where V c = N c /n is the roton cluster volume and the parameters o>f are given by Eq.fJHJ). 
Because Eq.(JIJ) follows from a Hartree-Fock variational procedure the number of atoms 
forming the roton cluster should satisfy the condition N c > 1 (not necessarily N c ^> 1). 
But below we will suppose the spherical symmetry of the roton clusters, which can only be 
accurate with a large number of atoms (for example one may demand N c > 10). 

Despite its apparently similar form, equation flU) should not be confused with the well- 

nn 

known Gross-Pitaevskii equation 8, 9|. Both contain the same hard-sphere repulsive poten- 
tial to describe short-range collisional interactions. However, in the Gross-Pitaevski equation 
the attractive harmonic potential, if present, describes an external trap; in equation (0} the 
harmonic Z7a(x) is a effective long-range many-body interaction, found by the variational 
procedure. Thus while the Gross-Pitaevski equation is only valid in the limit of low particle 
density, a^n 1, and small wavenumbers, ka>Q <C 1, equation (j3J) remains valid at high 
densities, a^n ~ 1, provided \k — fcol^o <^ 1 (o r more precisely, provided \k — ko\ < 0.2 A -1 ). 
In the limit when N c — > oo and uj\ — > the confining many-body potential {7a (x) vanishes 
and equation (HJ) does becomes just the Gross-Pitaevskii equation. 

Equation (jlj) is semiclassical, and as mentioned above, all phenomena in helium II, in- 
cluding roton excitations, should be described quantum mechanically. We present below a 
quantisation procedure for the NLSE based on a self-similar solution. 

Following the method of (for the 3D case see ^lj) it can be shown that this NLSE (j3J) 
has a self-similar solution when the dimensionless parameter e c = N^ 3 /[6(N C — 1)] satisfies 
the condition e c 1. Evidently, this is true for iV c > 10 because in this case e c < 0.04, 
and even for JV C = 3 we find that this parameter is small enough: e c = 0.12. Further, 



we have already shown |6( that the self-similar evolution of the atoms in the roton clusters 
described by the NLSE (j3J) can be derived from a Hamiltonian in the space of the ellipsoidal 
coordinates of the roton cluster: 

where n k = modat/dt are the canonical momenta, a k are the ellipsoidal coordinates, and 
mo = m/7 is the renormalized mass. The interaction constant G and the parameters io\ of 
the self-consistent confining potential C/r(oi, a 2 , a 3 ) = (m /2) J2k=i ^1°^ °f the roton cluster 
in Eq.(jSJ) are given by Q| : 



., 15(iV c - l)q ft 2 = 15(AT C - l ) ao n 
m^ia^a^a 2 , ' 7m 



The quantisation of the Hamiltonian (|SJ) yields a Schrodinger equation for the roton wave 
function xiflu a 2, a 3) '■ 

[H W = (-2^g^f + V(ai ' a2 ' G3) ) X ' (?) 
where the potential V(ax, 02,03) has the form 

7(01,02,03) = ^^^^+ ° ■ (8) 
z fc=1 010203 

The parameters uj\ must satisfy the stationarity conditions (dV/da k ) ak=3jk = 0, which just 
yield Eq.© again. We note that the equation (J7|) is not the standard Schrodinger equation 
because it given in the space of ellipsoidal parameters o&. This Schrodinger type equation 
can be solved by an expansion of the potential energy (jHJ) in a series with the variable 
Uk — a k — o-k, subject to the condition \yk\/cik 1- in the quadratic approximation we find 

7T&Q ^ ^ 

7(01,02,03) = v Q + — J2J2 A k S ykVs, (9) 

1 k=l s=l 

where the constant Vq and the matrix Ak s are given by 

5G 

Vo = n , A fcs = u k u s + 2u 2 k 5 ks . (10) 

20l0203 

The eigenvectors u^ k ' of the matrix A ks obey 



where the second equation is just the normalisation condition. Thus the coordinate trans- 
formation 
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y s = a s - a s = ^ uf ] z k} (12) 

k=l 

is orthogonal. Using Eqs.(P fT2"j) . we can now write the Schrodinger equation (jHj) as 

, t <9 / h 2 3 



+ (is) 

o k= i oz k l k=l ) 



dt \ 2m, 

We note that the energy spectrum of Eq. (fTBl does not depend on the renormalized mass m 
(which for some problems can be redefined ll|). Here the eigenfrequencies il^ are given by 
equation det(A — f2|I) = where I is the identity matrix, or, in explicit form, 

(pi - 3w?) (nj| - 3uj 2 2 ) (n 2 k - 3t4) - uolul (nj - 3^) 

-u\ul (p, 2 k - 3ufj - ujIu 2 2 (pi - 3ufj - 2ujIuj 2 2 ujI = . (14) 
This allows us to find the eigenfunctions Xn in2 n 3 and eigenenergies E nin2ri3 by 

(15) 

where the Hamiltonian 7i is that defined by the Schrodinger equation (J13)) . These eigen- 
functions Xnmans are 

Xmn 2 n 3 («l.«2.«3)= ^ ^ J ^ Wi!w 2 !w 3 ! J 1\ ™ ^~ 2™ ' ^ 



where H n (£) are Hermite polynomials and = ^ymo^fc/fo- The eigenenergies in Eq.(Jl 
have the form 

3 

S ninans = Vb + ^E Q fe( n fc + 1 /2), (17) 
fc=i 

where n& = 0, 1,2, ... at = 1,2,3 are the quantum numbers of a 3D quantum harmonic 
oscillator. Assuming now spherical symmetry of the roton clusters, we find at = a (for 
k — 1, 2, 3) where a is given by p] : 

a ^ I/3 - Hisf ■ < 18 > 

Here n = p/m is the particle density of liquid helium, giving q Q = 2.22 A. For example, in 
the case N c = 13 we have a = 5.22 A. We note that the value of the parameter go defined 
by the particle density n of liquid helium is very close to the value of the s-scattering length 



ao = 2.2 A for the helium atoms. From Eq. ffTHj) it also follows that a = go when N c = 1, so 
the physical sense of the parameter q is the same as that of a : an effective radius of the 
helium atoms. 

In this spherically symmetric case Eq. (jl4)) yields the eigenfrequencies f2i = O2 = y2uj, 
^3 = V^uJ. Here the frequency uj = [15(iV c — l)^] 1 ^ 2 ^?^ -1 ^ -5 / 2 follows from Eq.([H]l when 
afc = a for (k = 1, 2, 3). Let us define the energy of a symmetric excitation A s = hQ 3 and the 
energy of an asymmetric excitation A a = hfli = M7 2 - Then, taking into account Eq. ffTHj) . 
we find 



A. = 5^f^ {Nc ~,)r ^ A a = ^A s . (19) 



3 ROTON GAP 



At first we consider in detail the calculation of the s-scattering length ao for helium 
atoms based on the Lennard- Jones intermolecular potential. The Schrodinger equation for 
the s-scattering wave and appropriate boundary conditions are : 

0/77 

0"(r) + A; 2 0(r) = ^U{r)<j>{r), (20) 

0(0) = 0, 0(r) -> ^sin(A;r + 5 (k)) as r -> 00, (21) 

where = is the wavenumber, 0(r) is the radial s-scattering wavefunction, So(k) is 
the phase shift of the s-scattering wave, m c is the reduced mass, and U(r) is the scattering 
potential. Defining the s-wave scattering length as usual by Oq 1 = — lim fc ^ [/ccotg S (k)], 



we can write it in the form 



12] 



a = lim — j~ / U(r)(p(r)rdr, (22) 
fe^o % Jo 

(with the appropriate normalisation condition for the wave function). When the Born 
approximation is valid one can make in Eq. (J22|) the substitution 0(r) — > 0o( r ) where 
0o( r ) — k~ l sm(kr) is the wave function of free particle. However, for slow particles this 
approximation is generally meaningless because of divergence of the integral in Eq. (|2*2*|) . 
In the case ao > we can instead formulate a modified Born approximation (MBA), based 
on the substitution <p(r) — > 9(r — ao)fc _1 sin(Ar) in Eq. (j22|) . where 9{r) is the Heaviside step 
function. Thus we suppose that the wave function is zero when r < a (i.e. the region r < a 
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is unattainable for a slow particle). This leads to the equation : 

2 rfi roc 

a = / U(r)r 2 dr. (23) 

For example, let us consider intermolecular interactions with the potential U(r) = ar~ n 
where a > 0, n > 3. Eq.(j23jl then yields 

/2 mQ \ 1 /M) / 1 N i/(n-2) 



The exact solution of Eqs. (j2QI21j) for this potential leads to an equation which differs from 
Eq.fl2H) only in the coefficient: C n = (n-2)- 2 ^ n -^T(^)/T(^), where T(x) is the gamma 
function. Note that in the Born approximation the integral diverges in Eq.([22p. For example, 
at n = 6 we find C§ = 0.76 in the MBA and Cq = 0.68 in the exact solution; at n = 12 we 
have C\2 = 0.79 in the MBA and Cyi = 0.71 in the exact solution. 

We note that the MBA leads to the most accurate results when the potential U(r) 
has some characteristic length tq where ao < Tq. It asymptotically leads to the exact 
solutions for exactly integrable model potentials when a < r . The condition a < r does 
hold for helium intermolecular interactions given by the Lennard- Jones potential: U(r) = 
4e[(ro/r) 12 — (ro/r) 6 ]. In this case calculating the integral in Eq. (|23jl one can find a fifth 
order algebraic equation : 

p 5 -3p 2 -A = 0, \ = ^— 2 , (25) 

where p = (ro/ao) 2 and m c = m/2. We use in this paper the parameters of the intermolecular 
Lennard- Jones potential for helium atoms calculated by a self-consistent-field Hartree-Fock 
method e/k B = 10.6 K and r m = 2.98 A, where r m = 2^r . For these values we find 
from Eq.([25j) that p = 1.468, and a = 2.19 A. Taking into account the limited accuracy of 
our calculations we will assume in this paper that ao = 2.2 A jfl^l, which is very close to 
the result based on this MBA method. 

Let us denote £r the energy per particle of the roton cluster in the ground state and 
Eq the energy per particle of ground state of the bulk. Then taking into account that the 
energy of a bound state of two particles in the roton cluster is given by e = —U(r m ) (to 
within some additive constant e&) and the energy of two particles in the ground state of the 
bulk is e = —U (2a ) (to within an additive constant e&) we can find 

£r = -e - e 6 , S = -e - e b , (26) 



where e& is the contribution of the bulk to the ground state of the roton cluster. Here U(r) 
is the Lennard- Jones potential. Thus the chemical potential fi and the gap A of the roton 
cluster are : 

/j, = £ R -£ = -e + e , A = e - e , (27) 

The equation A = e — eo is fulfilled because from Eq.((TJ) it follows that ji = —A. Using 
the above definitions one can find the positive parameter eo = — 4e[(r /2ao) 12 — (r /2ao) 6 ] 
which yields e /k B = 1.95 K. Thus from Eq.(|2Tj) follows that A/k B = 8.65 K. This 
value coincides to three decimal places with the experimental value of the roton gap ?1: 
A/k B = (8.65 ± 0.04) K. Because £ Q is the energy per particle of the ground state of the 
bulk we can write 

E 

lim = £ = -e - e b , (28) 

where Eq is the ground energy of the bulk and N is the number of atoms in the bulk. The 
ground state energy per particle £q in the helium bulk can be found via the diffusion Monte 
Carlo method (l6lll7|: £o/ks = —7.15 K. This yields the bulk parameter e^/ks = 5.2 K. 



4 ROTON CLUSTER EXCITATION NUMBERS 

The equation A s = — \x holds only for the symmetric roton cluster excitation [6j with 
vibrational quantum numbers n\ = n 2 = and n 3 = 1. In the general case, the roton cluster 
excitations with arbitrary vibrational quantum numbers satisfy the generalised equation for 
the number of atoms N a in the roton cluster: E nin2n3 — E 000 = — /i, where E niTl2riA > Eq Q0 . 
Hence, combining Eqs.([17|) and (}2Tj) we find the generalised equation for the number of 
atoms N a in the roton cluster in the form 



n a A a + n 3 A s = e - e , (29) 

where n a = n\ + n 2 (with ri\ + n 2 + n% > 0), A s = and A a = hQ\ = hQ 2 - Using 
Eq. (fTI?|) we can write Eq. (j29|) in explicit form 



e - e 



\QoJ N? /b U2/5n a + n 3 )E a 



(30) 

where again n a + n 3 > 0, and E a = h 2 /(mal). Eq. pUj) is equivalent to a fifth order algebraic 
equation for the variable iV c . We can formulate the stability conditions of the roton cluster 
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FIG. 1: The full curve is the dimensionless function A S (N C )/E a (left side in (JSUJ)). The broken 
line is the constant (e — €o)/E a (right side in ()3U|) at n a = 0, = 1), and the chain line is the 
constant i/5/2(e — eo)/E a (right side in (|3T))) at n a = 1, n% = 0). This figure shows that a stable 
(dA s /dN c < 0) solution of Ea. (|3Uj) at n a = 0,ra3 = 1 is N c = Nqi = 13.1 ~ 13, and there are no 
solutions of Ea.(|3U|) at n a = 1,^3 = 0. It also shows that stable (dA s /dN c < 0) roton clusters are 
possible only when the integral numbers N c satisfy N c > 3. 

as d(E nin2n3 — E Q00 )/dN c < and Ad > 2a, where we are considering here the number of 
atoms iV c in the roton cluster as a variable parameter. Combining these equations with 
Eq.((7|) and Eq. ffTHj) one can rewrite these two stability conditions in the form 

dA f 



dN r 



< 0, \ D > 2q N 1 J 3 , 



(31) 



where A S (N C ) is given by Eq. (jl9|) and Xd = 12.6T _1 / 2 A is the thermal de Broglie wavelength 
of a helium atom. The numerical solutions of Eq. ()30|) for the parameter N c = N nan:i depend- 
ing on quantum numbers n a = n\ + ri2 and for a few lower-lying roton cluster excitations 
satisfying the the stability condition dA s /dN c < are (numerical results in parentheses): 
N 01 = 13 (13.1); N 20 = 28 (28.3); N u = 63 (62.7); N 30 = 99 (99.3); iV 02 = 117 (116.6); N 21 = 
170 (170);iVi2 = 268 (267.8); N 03 = 397 (397). The cases n a = 0, n 3 = 1 and n a = l,n 3 = 
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are shown in Fig. 1. 

The second condition in flHTJ) yields T < 1.5 K for N c = N 01 = 13; T < 0.9 K for 
N c = A^2o = 28; and, for example, T < 0.3 K for N c = iVo2 = 117. But it is known that the 
spectrum of phonon excitations in liquid helium II explains the Debye law of the heat 
capacity ~ T 3 when T < 0.6 K. Assuming that the symmetric vibrational excitations of the 
roton clusters are more probable then asymmetric we find that the stable roton clusters in 
the region 0.6 < T < 1.5 K in the main consist of 13 helium atoms. Eq. (j29j) at n a = and 
n 3 = 1 yields A s = A and hence N c = 13. 



5 PARAMETERS OF THE ROTON EXCITATION SPECTRUM 

One can calculate the value A s /ks = 8.67 K at N c = 13 which coincides within ~ 0.2% 
with our theoretical value A/ks = 8.65 K. The effective roton mass in Eq.Q is [6] fiR = 
pl/(4N c (JC)), where (K.) = (/Cj) + (/Ct) is the average energy of a particle in the roton 
cluster. Here the internal energy has the value p (JCi)/kB = 4.95 K at T = 1.1 K and 
(JC T )/E a = 0.96(T/T C ) 5 / 2 is the dimensionless energy of an ideal Bose gas per particle where 
E a = h 2 /{mafy. Employing in this last formula the experimental value T c = 2.18 K gives 
(K T )/k B = 0.43 K at T = 1.1 K. This yields ^ R /m = 0.16 for N c = 13 and k = p /h = 
1.92 A -1 , coinciding with the experimental value [7|. 

The momentum po at the minimum of roton spectrum in Eq.(^) can be approximately 



evaluated 



Q 



aspo — 2tt%/TZ^ where 7Z is the average distance between atoms in the localized 
bound state. In the classical approximation one can assume that 1Z = r m , so that k = 
Po/h = 2.1 A -1 . In the quantum case, for symmetric vibrational excitations (n\ = ri2 = 



0,n 3 = 1), Eq. ffTTj) yields £ 3 = z 3 /TZ q where 1Z q = y h / (m Q 3 ) , and from Eq. (j29|) we have 
that A s = A. Hence the momentum po at the minimum of the roton spectrum excitation 
has the form : 

2n% 



Po = = 27r^/m A, (32) 
where mn = m/7 and A = e — en. From this one can calculate kn = pn/h = 2 A -1 , about 

[in 

4% from the experimental value [3, One also can find a boundary point k c = p c /h 
in the elementary excitation spectrum between linear (phonon) and nonlinear regions as 
k c = 2ir/d c , where d c = 2a = 10.44 A is the diameter of the roton cluster at N c = 13. Hence 
k c = 0.6 A" 1 , which again agrees with experiment 0, Q| • We also note that the formfactor 
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S(k) for neutron scattering in liquid helium II is connected with the energy spectrum of 
roton cluster excitations e(p) (see Eq.flTJ) by Feynman's well-known approximate formula 
Q: S(k) = h 2 k 2 /{2me{hk)). Here 1.7A~ l < k < 2.1A~ l because Eq.© is valid only in 
the vicinity of the minimum of the roton spectrum excitation k = p /h. 



6 CONCLUSION 



In conclusion, we have proposed a full theoretical description of the model of roton 
cluster excitations in liquid helium II. Our work shows the important role of the attractive 
forces in the dynamic behavior of the roton excitations in helium that leads for some range 
of temperatures to the formation of localised bound states or clusters which also can be 
treated as quantum 3D solitons in the liquid helium II. It is shown that for some range 
of temperatures the smallest and most stable roton clusters consist of 13 helium atoms in 
symmetric vibrational quantum state. A natural model of this roton cluster has a central 
atom surrounded by a shell of 12 atoms situated at the vertices of a regular icosahedron ; 
the stability of this configuration is favored by its having the greatest number (six) of nearest 
neighbors for each atom in a shell. We have found theoretically all parameters defining 
Landau's roton excitation spectrum, in agreement to high accuracy with experimental data. 
For example, the roton gap given by the formula A = e — eo coincides to three significant 
figures (the experimental accuracy) with the experimental value measured by Yarnell et al []| 
from long-wavelength neutron scattering. 
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